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RESUMEN 

En estc trabajo modclamos la emision csperada de la lmca molecular 
C 17 0(J=3 — >2) en discos protoplanetarios, modificando diferentes parametros 
ffsicos para obtener distintas caracteristicas observacionales. Nuestra meta es de- 
tcrminar la clase de observaciones que nos permitiran extraer information sobre 
los parametros ffsicos de los discos. Con este proposito realizamos un analisis es- 
tadi'stico de componentes principales y una correlation lineal multiple en el conjunto 
de resultados obtenidos a partir de los modelos. Ademas, presentamos un estudio 
sobre futuras observaciones de lfnea molecular en discos protoplanetarios usando 
SMA y ALMA. 

ABSTRACT 

In this work we model the expected emission from the molecular line 
C 17 0(J=3 — >2) in protoplanetary disks, modifying different physical parameters 
to obtain distinctive observational signatures. Our aim is to determine the kind 
of observations that will allow us to extract information about the physical pa- 
rameters of disks. With this purpose we perform a statistical analysis of principal 
components and a multiple linear correlation on our set of results from the models. 
We also present prospects for future molecular line observations of protoplanetary 
disks using SMA and ALMA. 

Key Words: ISM: MOLECULES METHOD: STATISTICAL RA- 
DIATIVE TRANSFER STARS: PLANETARY SYS- 
TEMS: PROTOPLANETARY DISKS 



1. INTRODUCTION 
The mass of disks around young stars is ~ 99% 
gas and only 1% dust. However, since the dust opac- 
ity is large compared to the gas opacity in a wider 
range of wavelengths, the dust component domi- 
nates the absorption and reprocessing of stellar radi- 
ation and the emergent spectral energy distributions 
(SEDs) of disks around Classical T Tauri (CTTS) 
and Herbig Ae stars (HAe) . Most of the models con- 
structed to explain observed SEDs have taken into 
account only the dust component in the calculation 
of the disk temperature (Chiang & Goldrcich 1997; 
D'Alessio ct al. 1998, 1999, 2001). These models 
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have been successful to explain the observed SEDs. 
However, SEDs are not sensitive to kinematical infor- 
mation, details of the radial and vertical temperature 
distribution, the chemistry of the gaseous component 
(Aikawa ct al. 1996, 1997; Aikawa & Herbst 1999; van 
Zadelhoff ct al. 2001; Aikawa & Nomura 2006), the 
possibility that the gas in the upper layers is hotter 
than the dust in lower layers by absorption of UV 
and X-rays (i.e., Glassgold ct al. 2004; Jonkheid ct 
al. 2004; Kamp & Dullemond 2004) and/or photo- 
electrical effect (i.e., Nomura & Millar 2005), etc. 

The study of molecular line emission from disks 
around young stars is an important tool to infer 
physical characteristics of disks (e.g., Dartois et al. 
2003; Carr ct al. 2004; Pictu et al. 2005; Qi et al. 
2006; Raman ct al. 2006; Dutrcy et al. 2007). An 
advantage of a spectral line is that emission at differ- 
ent frequencies/velocities might be probing different 
disk regions, making lines an important test for disk 
models. However, this is also a disadvantage in some 
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sense, since one observes intensities convolved with 
the beam of the telescope and with a finite spectral 
resolution, and therefore the information of the dif- 
ferent regions is mixed up in a complex way. Thus, 
the analysis and relationship between observations 
and model properties might become very compli- 
cated and difficult to disentangle. There are some 
previous works that compare molecular line emis- 
sion from protoplanetary or circumbinary disks with 
specific models of such an emission (e.g., Koerner et 
al. 1993; Guilloteau & Dutrey 1998; Qi et al. 2003, 
2004) and, in general, agreement between model and 
observations is fairly good, at least in the general 
appearance of the maps. However, given the great 
deal of physical parameters involved in the resulting 
molecular line emission, it is not straightforward to 
determine those parameters from a particular obser- 
vation just by fitting an emission model. 

Our main aim in this paper is to identify a set of 
observational characteristics that give most informa- 
tion on the physical properties of the disk. Such ob- 
servational characteristics should then be given the 
heaviest weight in a fit between observations and 
model aiming to determine physical parameters in 
a disk. 

This paper is structured as follows: in §2 wc show 
the outline of this work. In §3 we describe the as- 
sumptions to calculate the disk structure models and 
the initial input parameters. In §4 we explain the 
radiative transfer calculation and we discuss the se- 
lection of the C 17 (J=3 — ► 2) as the molecular line 
to make our study. In §5 we outline the network of 
models and the general trends of the line emission 
maps. In §6 wc describe the statistical study that 
identify the set of observational characteristics that 
give most information on the physical properties of 
the disk, and we comment the results derived in §7. 
In §8 we show the comments and prospects for these 
studies. Finally, in §9 we perform a study of the 
detectability of our modeled disks with SMA and 
ALMA and we summarize the conclusions in §10. 

2. GENERAL CONSIDERATIONS 

In order to achieve the goal described in the in- 
troduction, we have developed a set of molecular 
line emission models calculated for various mass ac- 
cretion rates, radius, viscosities and maximum dust 
grain radius distributions. The small size scale of a 
protoplanetary disk (~100 AU) and their low tem- 
peratures (~ 100 K; Bcckwith et al. 1990; Miyake & 
Nakagawa 1995; Beckwith et al. 2000) require obser- 
vations with high sensitivity and subarsecond high- 
angular resolution, since 100 AU subtends 0'.'7 at 140 



pc (the distance to the Taurus cloud). This makes 
interferometric observations necessary to carry out 
this kind of studies. With the intention of repro- 
ducing a real interferometric observation of a proto- 
planetary disk with different physical parameters, we 
integrated the radiative transfer equation and con- 
volved each model with a beam of 0'.'4, as a compro- 
mise between resolution and sensitivity. From each 
resultant map, we have measured different observa- 
tional signatures, as if they were data from a real in- 
terferometric observation. Finally in order to obtain 
the best combinations of such observational param- 
eters that yield more information about the physical 
characteristics of disks, we have undertook a novel 
statistical approach to link observational properties 
of the expected molecular line emission with the un- 
derlying physical properties of the disk, by means of 
a principal component and multiple linear correla- 
tion analysis. We show that this is a promising type 
of analysis to prepare the observations with the new 
generation of millimeter and submillimcter interfer- 
ometers. 

In this study it is important to choose an appro- 
priated molecular transition sensitive to the physical 
parameters in which we may be interested. To carry 
out our study we have selected the C 17 0(J=3 — > 2) 
transition at 337 GHz. This line is a high excitation 
transition of a CO isotope with very low abundance, 
which makes it less susceptible to be affected by ab- 
sorption and/or the emission from the surrounding 
cloud material. This transition is also a suitable can- 
didate to be observed in protoplanetary disks us- 
ing the Submillimcter Array (SMA) and the Ata- 
cama Large Millimeter Array (ALMA) , as shown by 
Gomez & D'Alessio (2000). 

3. DISK STRUCTURE MODELS 

3.1. Assumptions 

We base our calculations of molecular line emis- 
sion on structure models of accretion disks irradiated 
by the central star, which have been previously used 
to explain different observations of classical T Tauri 
Stars. The assumptions and calculation method of 
such models are described in D'Alessio et al. (1998, 
1999, 2001). In summary, the disk is assumed to be 
in steady state, with a constant mass accretion rate 
M and an a— viscosity (Shakura & Sunyacv 1973), 
with a constant value of the viscosity parameter a. 
The disk is in vertical hydrostatic equilibrium in the 
gravitational potential well of the star, and we ne- 
glect the disk self-gravity. We assume that gas and 
dust are thermally coupled, having the same temper- 
ature everywhere. This dust/gas temperature enters 
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in the calculation of the disk volumetric density dis- 
tribution through the integration of the hydrostatic 
equilibrium equation. The main heating mechanisms 
considered are viscous dissipation and stellar irradi- 
ation. The viscous dissipation is important in heat- 
ing the inner regions (close to the star and close to 
the midplane) , the direct stellar irradiation heats the 
disk atmosphere, and the stellar radiation scattered 
and reprocessed by the disk upper layers heats the 
whole vertical structure. The transfer of radiation 
through the disk is calculated taking into account 
that the dust scatters and absorbs stellar and disk 
radiation, implying that the temperature structure 
depends on the dust properties. The viscous irradi- 
ated disk models used in the present study show the 
temperature inversion previously found by Calvet et 
al. (1991, 1992), i.e., at the outer disk, R > 10 AU, 
the upper layers are hotter than the disk midplane, 
because they are heated by direct stellar irradiation 
(see also Chiang & Goldrcich 1997 and D'Alessio et 
al. 1998). 

The dust opacity is calculated using the Mie the- 
ory for compact spherical grains. We consider a dis- 
tribution of sizes given by n(a) = noa~ p , where 
a is the radius of the grains, hq is a normaliza- 
tion constant, and p is a free parameter. In this 
work, we have adopted p = 3.5 from Draine & Lee 
(1984) and the model of dust composition proposed 
by Pollack et al. (1994) with the variations intro- 
duced by D'Alessio et al. (2001). To account for the 
possibility of dust growth, we adopt different val- 
ues of maximum grain sizes. Dust grains of different 
sizes have different continuum opacity at mm wave- 
lengths, affecting the molecular line emission in dif- 
ferent ways. The existence of bigger grains in disks 
than in the interstellar medium was proposed to ex- 
plain the observed slope of the continuum SED at 
millimeter wavelengths (Beckwith & Sargent 1991; 
Miyake & Nakagawa 1995). Large grains could be 
depleted from higher layers of the disk, but could be 
well mixed with gas below a few gas scale heights. 
For simplicity, the disk models adopted here (from 
D'Alessio et al. 2001), assume that gas and dust are 
well mixed. This seems a reasonable assumption if 
the emission from the molecular line arises mainly 
from areas closer to the midplane than to the upper 
layers. This point will be discussed later (subsec- 
tion 4.1). 

3.2. Input parameters 

For the present study we have considered the 
following input parameters: maximum disk radius 
(Rd), maximum radius of dust grains (a max ), disk 



mass accretion rate (M), and viscosity parameter 
(a) (see values in Table 1). 

TABLE 1 

INITIAL PHYSICAL PARAMETERS 





@>max 


M c 




(AU) 


(H 


(M©/year) 




50 


i 


io- 9 


0.001 


100 


10 


3 10~ 8 


0.005 


150 


10 2 


io- 7 


0.01 




10 3 




0.02 




10 4 




0.05 




10 5 







a Disk radius. 

b Maximum radius of dust grains. 
c Mass accretion rate. 
d Viscosity parameter. 



We have adopted typical parameters of a T Tauri 
star from Gullbring et al. (1998), i.e., M* = 0.5 M , 
i?* = 2 Rq, and T* = 4 000 K for all the models. The 
disks are assumed to be at 140 pc, the distance of the 
Taurus molecular cloud (Kenyon et al. 1994), with a 
typical inclination angle i = 60°. It is important to 
mention that each disk structure is sclf-consistently 
calculated given these input parameters. This means 
that the whole disk structure is affected by all the pa- 
rameters, consequently affecting the line properties. 
This might complicate the analysis of the resulting 
line properties, but we think this gives a more realis- 
tic description of the interplay between the different 
variables. 

4. RADIATIVE TRANSFER AND MOLECULAR 
LINE EMISSION 

The model of the disk structure provides a de- 
tailed density and temperature distribution through 
the disk as a function of the height and the distance 
to the disk center. To derive the line intensity for 
a given molecular transition, we must resolve the 
transfer equation. We have used the same assump- 
tions and formalism that Gomez & D'Alessio (2000), 
and we summarize them briefly here. We assume 
local thermal equilibrium for the population of the 
molecular energy levels and we consider thermal line 
profiles. We divide the disk in a grid of cells consid- 
ering isovelocity lines and their perpendicular lines 
(see Appendix in Gomez & D'Alessio 2000). We in- 
tegrate the transfer equation dl v /ds = K v p{S l/ — I u ) 
through the line of sight at the center of each cell, 
where /„ is the intensity, s is the length along the line 
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of sight, k u is the absorption coefficient that consider 
the contributions from the line and continuum, i.e., 
K v = ki + k c , p is the mass density of the gas, and S v 
the source function. The coefficient k c is dominated 
by dust and we consider pure absorption opacity of 
each kind of dust grains size. 

In order to derive the flux density, we convolved 
all our models with a beam of 0'.'4 of HPBW, as a 
compromise between resolution and sensitivity. Fi- 
nally we subtract the continuum emission to isolate 
the flux density of the molecular line, obtaining a 
set of model results that reproduce a real observa- 
tion of a protoplanctary disk with different physical 
parameters. 

4.1. Selection of the emission line transition 

It is very common that young stellar objects like 
T Tauri stars are still embedded in the material of 
the parental cloud. The envelope that surrounds the 
disk-star system is composed by cold gas and dust 
that could hide the disk emission. 

On one hand, the surrounding material could ab- 
sorb the emission from the hotter, innermost part 
of the disk structure. This problem can be solved 
by selecting a molecular transition whose energy is 
high enough to trace the hotter gas from the disk, 
while few molecules in the colder envelope are in the 
states involved in this transition, thus avoiding line 
absorption at the envelope. However, the frequency 
of the transition must not be too high, because oth- 
erwise the dust in the envelope would become opti- 
cally thick. On the other hand, the molecular line 
emission from the whole cloud could hide the deeper 
emission from the disk, as it is usually the case for the 
more abundant CO isotopes. To make sure that the 
observed emission comes only from the disk struc- 
ture, we must select a molecular species of low abun- 
dance for which the envelope is optically thin and 
therefore its emission would be negligible. 

In order to fulfill these conditions, we have chosen 
the (J=3 — > 2) transition of the C 17 molecule at 
337 GHz. Its molecular abundance relative to H2 in 
molecular cores is low, 5.0xl0 -8 (White & Sandell 
1995) and the frequency transition at 337 GHz is 
still low enough to avoid being much absorbed by 
the dusty envelope. 

In a recent work, Dartois ct al. (2003) have stud- 
ied the disk vertical temperature structure using dif- 
ferent isotopes of CO. Being characterized by differ- 
ent opacities, different lines trace the gas at distinct 
depths. Dartois et al. (2003) find a good agreement 
between the inferred temperature for each transi- 
tion/isotope and the temperature where r ~ 1 in 



irradiated (dusty) disk models. Following Dartois et 
al. (2003), we have studied the formation regions of 
molecular lines from different isotopes of CO (CO, 
13 CO, C 18 0, and C 17 0) in our set of disk models, 
resulting that the C 17 is formed closer to the mid- 
plane than the rest of the isotopes. This fact makes 
the LTE assumption acceptable, since the region 
that mainly contribute to the line emission shows 
higher densities than the critical density of the se- 
lected transition (~ 5 x 10 4 cm -3 ). 

In addition, given the high gas densities in the 
disk, we expect that the abundance of this isotope, 
would be less affected by photodisociation produced 
by the incident radiation, since it would be shielded 
against it, specially in the deepest layers. More- 
over, theoretical models of the evolution of molec- 
ular abundances in protoplanctary disks, predict de- 
pletion of CO from the gas phase for temperatures 
below 20 K (Aikawa et al. 1996) and therefore, de- 
pletion is probably not significant for the chosen 
molecule and for the disk radii we are considering 
here (< 150 AU). Thus, for simplicity, we adopt a 
constant abundance for C 17 relative to H2, given 
the molecular core value (White & Sandell 1995). 



5. RESULTS 

We have calculated the expected emission in the 
C 17 0(J=3 — > 2) transition for a typical T Tauri star 
with a disk inclination angle of 60°, and for all pos- 
sible combinations of the physical parameters shown 
on Table 1. We have 241 different models for which 
we have integrated the radiative transfer equation 
at 12 different velocities, from -2.5 to 2.5 km s _1 
at steps of 0.5 km s _1 . The results are maps like 
those shown in Fig 1, in which we only show the 
positive velocities, since the maps are similar and al- 
most symmetrical in negative velocities to the ones 
shown here (there are slight differences though, due 
to the hypcrfine structure of C O transition). From 
each map we have measured the following observa- 
tional signatures as if they were data from a real 
interferometric observation: intensity of the princi- 
pal (the more distant from the observer, to the north 
in our maps) and secondary (closer to the observer, 
to the south in our maps) peak at each velocity, dis- 
tance from the disk center to principal peaks, half 
power size of the emission distribution, and velocity 
at which the maximum intensity is present. These 
represent a total of 43 different observational param- 
eters for each input model. 
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Fig. 1. Emission maps at different velocities for a disk 
with radius = 150 AU, a mal = 10 fim, i = 60°, a = 0.01 
and mass accretion rate = 1O _7 M0 year" 1 . Maps have 
been convolved with a 0"4 beam. The lowest contour 
and the increment step are 20 mjy beam -1 . The ellipse 
traces the outer edge of the disk. 

5.1. Line emission maps 

In all maps obtained from our set of models we 
observe the same tendencies as Gomez & D'Alessio 
(2000). Summarizing, we observe an asymmetry on 
both sides of the major axis, as expected from op- 
tically thick emission. The areas further away from 
the observer (positive declination in Fig. 1) show line 
emission of higher intensity because the line of sight 
intercepts areas of the disk closer to the central star, 
where the gas is warmer. This result is confirmed in 
DM Tau spectral line observations by Dartois et al. 
(2003), where an inclination angle i ~ —37° was as- 
sumed. We notice another asymmetry on both sides 
of the minor axis, that it is more pronounced when 
approaching to the velocity of the cloud (v = km 
s -1 ). It is caused by the asymmetry of the hyperfine 
structure in the C 17 molecular transitions. 

We find that the maximum intensity of the line at 
the systemic velocity (v = km s -1 ) traces the outer 
edge of the disk. As previously discussed by Sargent 
& Beckwith (1991) and Gomez & DAlessio (2000), 
this is a consequence of the fact that the effective 
area emitting at a given velocity within the beam 
increases with the distance to the central star more 
steeply than the decrease of brightness temperature 
with distance. 

Moreover, at the center of the disk the emission 
intensity diminishes, and we can even sec absorption 
lines in some of our maps. The high dust continuum 
opacities typical in the central parts of the disk, re- 
duces the contrast between the emission lines and 
the continuum, and, if it is high enough, lines could 
show up in absorption. 



6. STATISTICAL TOOLS 

As we mentioned in section 5, we chose a set of 
43 different observational parameters to character- 
ize the maps resulting from our model calculations. 
These parameters are, in principle, somewhat arbi- 
trary. If we want to simulate a real observation, in 
which we would like to extract information about the 
underlying physical characteristics of disks, it is obvi- 
ous that some of these 43 parameters will have more 
informative power, while some may turn out to be ir- 
relevant. Moreover, it is likely that not all the chosen 
observational parameters will be independent. 

Here, we have undertook a statistical analysis to 
try to identify a set or a combination of observa- 
tional parameters that could render more informa- 
tion about the disks properties. First, we used a 
principal component analysis to reduce the number 
of observational parameters to a small, informative 
set. Later, we investigate whether we can obtain 
quantitative values of physical magnitudes from ob- 
servational parameters, by means of multiple linear 
regression. 

6.1. Principal Components 

The principal component analysis is a statistical 
technique that provides a dimensional reduction of 
a set of variables. In our case, each of the initial 43 
observational parameters would be an axis in a sys- 
tem of coordinates in a multidimensional space. The 
method consists of finding a set of orthogonal axes 
in which the variance (heterogeneity) of our data is 
maximum (see Thurstone 1947; Kaiser 1958). This 
is solved through a linear and orthogonal transfor- 
mation that corresponds to a rigid rotation of the 
original data into a new set of coordinates. The 
eigenvalues obtained provide information about the 
variances of the data in the new space, and the eigen- 
vectors represent the direction of the axes in the new 
space of representation of our data. In our study, this 
analysis reduces the number of observational signa- 
tures necessary to derive information about the phys- 
ical parameters, to a smaller set of linearly indepen- 
dent parameters: the principal components. 

6.2. Multiple linear correlation 

In order to quantitatively estimate each physical 
parameter from the set of observational variables, we 
have also carried out a multiple regression analysis 
(Pearson 1908). In our case, the obtained principal 
components will be considered the independent vari- 
ables, and the physical properties of the disk will be 
the dependent ones. The multiple linear correlation 
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analysis will then try to make the best fit, to derive 
a linear function of the form: 

Ui = a\X\ + a 2 x 2 + ••• + ajXj , 

where yi is the dependent variable, and the 
independent ones. In an ideal case (correlation coef- 
ficient ~ 1), we would obtain a function with a good 
predictive power for the dependent variable. For our 
particular problem, we would like to obtain a formula 
with which, from a set of observational parameters, 
we could calculate the physical characteristics of the 
disk. 

Moreover, to test the signification of the multiple 
regression we applied the F-Snedecor test (Snedecor 
1934) as hypothesis testing. The ratio of two chi- 
squares divided by their respective degrees of free- 
dom follows an F distribution. The test consist of 
comparing the relation between the variance of the 
predicted values for the dependent valuable and the 
error variance with the value of the F distribution. 

7. OBSERVATIONAL SIGNATURES VS. 
PHYSICAL PARAMETERS 

7.1. Eigenvalues and Eigenvectors of Principal 
Components 

To reduce the number of relevant components in 
the principal component analysis, we adopted the 
Kaiser criterion (Kaiser 1960), which only retains 
factors with eigenvalues greater than 1. In our case, 
the criterion selects four factors, which account for 
91% of the total variance of the system. In fact, the 
first and the second factors alone represent 82% of 
the total variance. Therefore, most of the following 
analysis is based on these two first principal compo- 
nents. 

We have represented the eigenvectors in Table 2, 
for each principal component. The numerical entries 
in this table show the linear combination coefficients 
for each observational parameter, used to build the 
corresponding principal component. Therefore, they 
indicate the relative weight of each observational pa- 
rameter on the components. 

For the first principal component (PCI) the ob- 
servational parameters with a larger weight in its def- 
inition are (in order of decreasing relative weights) 
the velocity of the peak emission, the half power sizes 
for principal peaks at intermediates velocities, and 
the distance from principal peaks to center. The pa- 
rameters that define the second principal component 
(PC2) are the half power sizes of principal peaks at 
1.5 km s _1 , the velocity of the peak emission, and 
the half power sizes of secondary peaks at 0.0 and 
1.5 km s . The third principal component (PC3) is 



TABLE 2 
EIGENVECTORS 



Observational 
Parameters 


PCI 


PC2 


PC3 


PC4 


I(vo) a 


-0.0279 


0.0020 


0.0030 


0.0033 


I(V0. 5 ) a 


-0.0276 


0.0019 


0.0029 


0.0034 


I(vi) a 


-0.0281 


0.0007 


0.0022 


0.0035 


I(vi.5) a 


-0.0340 


-0.0023 


0.0024 


0.0040 


I(v 2 ) a 


-0.0276 


-0.0086 


0.0016 


0.0040 


I(v 2 .5) a 


-0.0118 


-0.0082 


0.0028 


0.0040 


r(vo) b 


-0.1360 


0.0388 


-0.0072 


-0.0137 


r(v . 5 ) b 


-0.1462 


0.0325 


-0.0041 


-0.0169 


r(vi) b 


-0.1581* 


0.0448 


-0.0015 


-0.0086 


r(vi. B ) b 


-0.1930 


0.0100 


0.0266 


-0.0262 


r(v 2 ) b 


-0.1077 


0.0075 


-0.0030 


-0.0150 


r(v 2 . 5 ) b 


-0.0483 


0.0310 


-0.0094 


-0.0062 


Isec(v )° 


-0.0289 


-0.0047 


0.0004 


0.0019 


Iaec(v .5) C 


-0.0297 


-0.0042 


0.0006 


0.0022 


WV!) C 


-0.0179 


-0.0091 


-0.0024 


-0.0018 


I«bc(vi.b) C 


-0.0003 


0.0007 


-0.0009 


-0.0004 


Wv 2 ) c 


0.0000 


0.0000 


0.0000 


0.0000 


Iaec(v 2 .5)° 


0.0000 


0.0000 


0.0000 


0.0000 


a(v ) d 


-0.0595 


0.0426 


0.0316 


0.0034 


b(v ) d 


-0.0790 


0.0334 


0.0272 


-0.0184 


a(v . 5 ) d 


-0.1615 


-0.0309 


0.0259 


0.0043 


b(v . 5 ) d 


-0.1947 


-0.0788 


0.0523 


-0.0076 


a(vi) d 


-0.1138 


-0.0250 


0.0281 


-0.0024 


b(vi) d 


-0.0388 


-0.0024 


0.0041 


0.0017 


a(vi. 5 ) d 


-0.1777 


-0.1926 


-0.0350 


-0.0615 


b(vi. 5 ) d 


-0.2271 


-0.1499 


0.0086 


0.0743 


a(v 2 ) d 


-0.2420 


-0.0171 


0.0543 


-0.0023 


b(v 2 ) d 


-0.1181 


0.0056 


0.0133 


-0.0084 


a(v 2 . 5 ) d 


-0.0286 


0.0168 


-0.0054 


0.0011 


b(v 2 . 5 ) d 


-0.0028 


0.0114 


-0.0067 


0.0010 


a scc (vo)° 


-0.0746 


0.1023 


0.0703 


0.0434 


Usee V v O j 


-0 0436 


0. 1026 




-O 0607 


a S ec(v .5)° 


-0.0026 


0.0161 


-0.0127 


-0.0076 


b se e(Vo. 5 ) 


-0.0006 


0.0047 


-0.0057 


-0.0011 


a sec (vi)° 


0.0000 


0.0000 


0.0000 


0.0000 


b sec (vi) e 


0.0000 


0.0000 


0.0000 


0.0000 


a»ec(vi.s) e 


-0.0733 


0.0985 


0.0683 


0.0525 


h„ c (vi.5) e 


-0.0398 


0.0951 


0.0398 


-0.0549 


a sec (v 2 ) c 


-0.0040 


0.0185 


-0.0141 


-0.0097 


b 3C e(v 2 )° 


-0.0008 


0.0043 


-0.0060 


-0.0006 


a SC c(v 2 . 5 )° 


0.0000 


0.0000 


0.0000 


0.0000 


b SC c(v 2 . 5 ) e 


0.0000 


0.0000 


0.0000 


0.0000 




0.3176 


-0.1431 


0.1527 


-0.0241 



1, 1.5, 2, 2.5 km s _1 ) in Jy beam -1 . 



b Distance from disk center to principal peaks in arcsec. 
intensity of the secondary peak at each velocity. 
d Minor (a) and mayor (b) axis of the half power sizes 
of emission for the principal peak of intensity at each 
different velocity, in arcsec. 

c Minor (a) and mayor (b) axis of the half power sizes 
of emission for the secondary peak of intensity at each 
different velocity. 

'Velocity at which the maximum intensity is present in 
km s~ . 

Values in boldface represent the observational parame- 
ters with a larger weight in the definition of each PC. 
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defined by the velocity of the peak emission and the 
half power sizes for principal peaks at 0.5 and 2.0 km 
s -1 , and for secondary peaks at 0.0 and 1.5 km s _1 . 
Finally, the physical parameters that define the forth 
principal component (PC4) are the half power sizes 
for principal peaks at 1.5 km s _1 and for secondary 
peaks at 0.0 and 1.5 km s _1 . We have marked in 
boldface these most representative observational pa- 
rameters in Table 2. 

7.2. PCI (kinematical component)-PC2(spatial 
component) diagrams 

Due to the fact that the first and the second prin- 
cipal component represent the 82% of the total vari- 
ance of the system, we have represented all our mod- 
els in a PC1-PC2 diagram, to check whether such 
diagrams can be used to discriminate disks with par- 
ticular physical characteristics, if we see some clus- 
tering related to the physical properties. Consider- 
ing the parameters with the larger weights in each of 
these principal components, we have named them as 
"kinematical component" and "spatial component" 
in the case of PCI and PC2, respectively. 

The most evident result is shown when 
we represent disks with different radii in the 
PCl(kinematical)-PC2(spatial) diagram (see Fig 2). 
The disks with radius 50 AU are distributed on the 
right part of the diagram, the ones with radius 100 
AU are located at the central part and the disks with 
radius 150 AU are located on the left of the plot. 
This means that the kinematic principal component 
(x axis) is qualitatively good to discriminate among 
disks with different radii. This result indicates that 
we can get information about the radius of a pro- 
toplanetary disk (physical parameter) from the ve- 
locity of the peak emission, the half power sizes for 
principal peaks at intermediates velocities and from 
the distance from principal peaks to the center (ob- 
servational signatures). That the disk radii can be 
discriminated relatively well by maps of line emis- 
sion may seem a relatively obvious result, but it il- 
lustrates the power of this kind of statistical analysis. 

Another interesting trend is seen from the repre- 
sentation in the diagram of disks with different mass 
accretion rate (see Fig 3). Disks with higher mass ac- 
cretion rates tend to populate the upper parts of the 
PCl(kincmatical)-PC2(spatial) diagram. Therefore, 
in this case it is the second principal component (y 
axis) the one that better discriminates among disks 
with different mass accretion rates. Considering the 
parameters that give rise to this component, we can 
say that the half power sizes of principal peaks at 
1.5 km s , the velocity of the peak emission and 



the half power sizes of secondary peaks at 0.0 and 
1.5 km s _1 velocities provide information about the 
mass accretion rate of the disks. 

Other trends relating principal components and 
the rest of the physical parameters are also present, 
but qualitatively they are not as clear as the two we 
have mentioned. In the case of other physical param- 
eters, it is the second principal component the one 
that provides more information about the a parame- 
ter and the maximum radius of dust grains. Further 
analysis including models of more molecular transi- 
tions will certainly be useful to obtain refined prin- 
cipal components that can better discriminate these 
physical parameters. 

7.3. Multiple correlation 

We have seen in the previous section that dia- 
grams of principal components can be useful to dis- 
criminate among disks with different physical char- 
acteristics. It would be interesting, however, to have 
a mathematical tool to easily obtain numerical val- 
ues for these physical characteristics based on the 
observed maps. As a first approach, we have here 
performed a multiple linear correlation with the ob- 
servational parameters as independent variables, and 
the physical characteristics as dependent ones. The 
resulting coefficients for such a fit are shown in Ta- 
ble 3. 

As it was deduced in the previous subsection, 
the linear correlation coefficients of each principal 
component (see r(PCj) values in Table 3) show 
that PCI (kinematical component) is the best one 
to derive information about the disk radii and 
PC2(spatial component) provide the most informa- 
tion in the determination of the rest of the physi- 
cal parameters, M , a, and maximum radius of dust 
grains. 

The most important piece of information that can 
be retrieved with this kind of study is the deduction 
of the linear combination of observational signatures 
that provide quantitative information about physi- 
cal characteristics. This linear combination can be 
derived from A^ coefficients in Table 3 as follows: 

P = A + AiPCl + A 2 PC2 + A 3 PC3 + A 4 PC4 

where P is the physical parameter. As an example, 
in the case of the radius of the protoplanetary disk, 
this combination would be the following: 

R d = -55.3 - (79.4 x PCI) - (54.2 x PC*2)+ 
+ (133.3 x PC3) - (137.7 x PC4) 
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TABLE 3 

CORRELATION COEFFICIENTS AND F-SNEDECOR TEST VALUES 







r(PCl) a 


r(PC2) a 


r(PC3) a 


r(PC4) a 


R b 


F 




Radius 


-0.9 


-0.19 


0.07 


-0.001 


0.97 


951.62 




M 


-0.18 


0.5 


-0.3 


-0.16 


0.57 


26.93 




a 


0.16 


-0.21 


0.13 


0.07 


0.31 


5.94 


Max. 


Grain Radius 


-0.08 


-0.009 


-0.012 


0.08 


0.19 


2.07 






A C 




A 2 C 


A 3 C 


A 4 C 






Radius 


-55.3 


-79.4 


-54.2 


133.3 


-137.7 






M 


1.8xl0~ 7 


-0.16xl0 -7 


1.4xl0~ 7 


-1.7xl0 -7 


-0.15xl0 -7 






a 


-0.0195 


0.0024 


-0.0349 


+0.0465 


-0.0091 




Max. 


Grain Radius 


-16943.4 


-100.4 


-40213.2 


+48703.1 


+123751.8 





a r(PCi) are the linear correlation coefficients of each principal component. 
b Multiple correlation coefficient. 

c Ai are the coefficients of the linear combination of principal components obtained from the multiple regres- 
sion study, to derive the values of the physical parameters in the first column. 



With these sets of linear combinations, we could 
estimate physical parameters from observations, pro- 
vided that the fit is good enough. 

The strongest correlation (see R values in Ta- 
ble 3) is obtained for radius (correlation coefficient R 
= 0.97), followed by mass accretion rate (R = 0.57), 
viscosity parameter (R = 0.31) and maximum radius 
of dust grains (R = 0.19). The multiple correlation 
coefficients alone are not good statistical indicators 
of the goodness of linear fits. To assess the validity 
of our fitted functions we carried out an F test to 
derive the signification of the linear regression. We 
assume a confident limit of 95%, which means that 
we could admit as good fits with values F < 1.44. 
However in all our results the values exceeded this 
critical value (see Table 3). This result suggests that 
our variables are far from the linear regime, which is 
certainly reasonable. 

7.4. Effects of measurements errors 

In our method, the derived value of a physical 
parameter depends linearly on the principal compo- 
nents, which are linear combinations of the observa- 
tional data. Therefore we must consider on one hand 
the error in the calculation of the principal compo- 
nents, and on the other hand the error derived from 
the calculation of the multiple correlation. 

In a real observation, the error of the intensity 
will depend on the rms noise reached during the ob- 
servations, and, in particular, the achieved signal- 
to-noise ratio will determine the positional accura- 
cies (ltr error in position ~ HPBW/(2 SNR)), as- 
suming a compact emitting source). In our method, 



for the determination of the most significant princi- 
pal components, the most important sources of er- 
ror are related to the observational parameters that 
contribute with a larger weight, i.e., the half power 
sizes of emission for the principal and secondary 
peaks of intensity (whose error can be roughly es- 
timated as ~ HPBW/SNR), and the velocity at 
which the maximum intensity is present (with la 
error ~ Av/ (2 SNR), where Av is the lincwidth). 
Since principal components are linear combinations 
of observational parameters, the weight of the er- 
rors associated with each measured observational pa- 
rameters will depend on theirs particular weight in 
the definition of each principal component (see Ta- 
ble 2). Nevertheless, we must point out that, in 
the case of calculations of the physical parameters 
of disks, the error is dominated by the uncertain- 
ties in the fits of the multiple regression. Since our 
derived variables are far from linear regime (see F 
test above) , measurements errors are not significant 
in comparison. 

8. COMMENTS AND PROSPECTS FOR THESE 
STATISTICAL STUDIES 

The results derived from the statistical study pre- 
sented in this paper show that this method could be 
a powerful tool to obtain information related to the 
physical characteristics in protoplanetary disks from 
observational parameters. 

The most important information of our study is 
related to the radius of the disks. This preliminary 
study results in a way to obtain a first approximation 
of the disks radii that depends only on the obser- 
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R=50 AU 



PCI (Kinematical component) 



-0.5 0.0 0.5 




-1.0 -0.5 0.0 

PCI (Kinematical component) 



R=100 AU 



PCI (Kinematical component) 



I.O -0.5 0.0 0.5 




-1.5 -1.0 -0.5 0.0 

PC 1 (Kinematical component) 



—i -0.4 



R=150 AU 



* : 



1.5 -1.0 -0.5 0.0 0.5 



PCI (Kinematical component) 




-1.0 -0.5 0.0 

PCI (Kinematical component) 



Fig. 2. Disk radii in the PC1-PC2 diagrams. Dots rep- 
resent all the modeled disks with different physical pa- 
rameters. Squares represent disks with radius 50 A.U. 
(upper), 100 A.U. (central) and 150 A.U. (bottom). 



Fig. 3. Mass accretion rates in the PC1-PC2 diagrams. 
Dots represent all the modeled disks with different phys- 
ical parameters. Squares represent disk with mass ac- 
cretion rate of 10~ 9 M^/year (upper), 3xl0 -8 Mo/year 
(central) and 10 -7 M^/year (bottom). 
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vations, avoiding the application of \ 2 fitting tech- 
niques. 

Since this is a promising method of study, and 
considering the future set of observations in proto- 
planetary disks that will be carried out with the de- 
velopment of ALMA, we plan to complete and im- 
prove our models and analysis. 

Nowadays it is not possible to test observation- 
ally our method because there is no observation of 
protoplanetary disks in our selected molecular tran- 
sition, C 17 0(J=3— >2). But, considering the excel- 
lent characteristics of this transition/isotope to ob- 
serve protoplanetary disks (see subsection 4.1), one 
should be able to test our study in the near future 
with ALMA observations. Meanwhile this moment 
comes, we will improve or study. 

First, to better determine the tendencies in the 
principal component diagram, we plan to increase 
the sample, building a more extended set of models 
with a wider variety of initial physical and obser- 
vational parameters. Calculations in other molec- 
ular transitions or isotopes would further constrain 
the information on temperature and density. An up- 
date in the calculation of the disk structure models 
would be done, considering in detail the photodiso- 
ciation or depiction effects in the line of study, along 
with dust grain growth and settling. We also plan 
to study the application of our method towards non- 
axial-symmctric disks, for what the development of 
a new disk structure that takes into account the lack 
of symmetry is needed. Finally, the determination 
of a more complex method to derive the physical 
parameters of disks from the principal components 
is desirable, since the application of multiple corre- 
lation provides variables far form the linear regime 
and, introducing the greater source of error in our 
study. 

9. DETECTABILITY OF MOLECULAR LINES 
WITH NEW INTERFEROMETERS 

It would obviously be important to be able to test 
our results observationally in the future. With this 
mind, we have also studied the detectability of the 
calculated models, when observed with interferome- 
ters that will be able to reach subarcsecond resolu- 
tion at the frequency of the C 17 0(J=3 — > 2) transi- 
tion, i.e., SMA and ALMA. 

We have calculated the line intensity expected for 
all our line models. The lowest peak intensity (~7 
mJy beam -1 ) corresponds to disks with Rd=5Q AU, 
M = 10 -9 M© yr -1 , a = 0.05 and 1 (im maximum 
grain radius. On the other hand, the highest peak 
intensity (~180 mJy beam -1 ) corresponds to disks 



with 150 AU radius, M = 10 -7 M© yr -1 , a = 0.01 
and 10 5 /jm maximum grain radius. 

With the SMA, a sensitivity of ~85 mJy beam -1 
is expected 4 , considering the 8 antennas of the array 
working at 337 GHz, with 1 km s -1 velocity res- 
olution and 10 hours integration time, under stan- 
dard values of precipitable water vapor (~2.0 mm 
for 300-355 GHz). 

On the other hand, in the case of ALMA, consid- 
ering 64 antennas, 1 km s -1 velocity resolution, only 
1 hour integration time and 1.5 mm of precipitable 
water vapor (median at the site over all hours and 
seasons), the sensitivity expected is 1.2 mJy beam -1 
(Butler & Wooten 1999). 

Considering 5er emission as detections, we con- 
clude that all our disks modeled, even the faintest, 
will be detected with one hour integration time with 
ALMA. Nevertheless, it will be extremely difficult 
of observe any of our modeled disks with the SMA, 
at least the ones with the physical characteristics 
showed in this work (see Table 1), since the SNR ob- 
tained for the modeled disks with the highest inten- 
sities (~180 mJy beam -1 ) in 10 hours of integration 
time, would only be 2a. 

Therefore, ALMA will be a crucial instrument to 
observationally test the predictions and assumptions 
of our models. 

10. CONCLUSIONS 

In this paper we have presented a statistical 
method to derive physical parameters from obser- 
vational characteristics in a protoplanetary disk. To 
make this study we modeled the expected emission 
of the C 17 0(J=3 — > 2) transition from protoplan- 
etary disks with different physical properties. We 
then applied a principal component and a multiple 
correlation analysis, to obtain a set of linear com- 
binations of observational parameters that may pro- 
vide relevant information of the physics of the disks. 
The main conclusions are the following: 

• The most significant results of our analysis are 
related to disk sizes. We can discriminate 
among disks with different radii using a princi- 
pal component composed mainly of the velocity 
of the peak emission, the half power sizes for 
principal peaks at intermediates velocities and 
the distance from principal peaks to center. 

• Moreover, some information about the mass ac- 
cretion rate could be obtained from a principal 
component made of the half power sizes of prin- 
cipal peaks at 1.5 km s -1 , the velocity of the 

4 http: / /smal. sma.hawaii.edu/call. html 
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peak emission and the half power sizes of sec- 
ondary peaks at 0.0 and 1.5 km s _1 velocities, 
although the results are much less significant 
than in the case of radii. 

• Although with preliminary results, the statisti- 
cal method presented here seems to be promis- 
ing and useful, and will be improved and com- 
pleted in a future with studies in other transi- 
tions. 

• We have performed a study of detcctability with 
ALMA and SMA. All our modeled disks could 
be detected with ALMA, using one hour integra- 
tion time, nevertheless the sensitivity reached 
by the SMA is not enough to detect our disks 
with reasonable integration times. We conclude 
that ALMA will play an important role to test 
observationally our models and our statistical 
results. 
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CyT and PAPIIT, Mexico. The model calculations 
were performed in the cluster at CRyA-UNAM ac- 
quired with CONACYT grants 36571-E, 47366-F, 
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